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jZ^, Abstract. Solutions to the collinear three-body problem which do not end in 

triple collision pass through an infinite number of binary collisions. Given three 

^L^ ■ masses, we show that four geometric quantities generate a finite description of 

itineraries of binary collisions. In the best circumstances, this description is 
semi-conjugate to a Poincare map of the flow. For other cases these quantities 
give upper and lower bounds on the itineraries which can occur. In addition 

rf\ ' to describing the dynamics of the collinear three-body problem, the results of 

Qthis paper rederives the existence of oscillatory motion in the ./V-body problem 
for N > 3. 
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1. Introduction 



£N| ' In the seventeenth century, Newton formulated the universal law of gravitation 

and completely solved the two-body problem. Moreover, his methods confirmed 
Kepler's laws of planetary motion. After this grand success Newton turned his 
fluxions to the motions of the Earth, moon and sun. He eventually gave up working 



£\i ■ on this three-body problem, saying it gave him headaches 1 10 

f^ , In the intervening years much work was done on the three-body problem gen- 

erating new techniques and new questions. After Poincare showed that chaotic 
behavior can occur in the restricted three-body problem p3 , the search for a com- 
plete solution halted abruptly. Poincare's result changed the issue from trying to 
solve an initial value problem to asking what behaviors can occur in the three-body 
C^ problem. 

A critical step in understanding what dynamics can occur in the three-body 
problem was identifying behavior near triple-collision. McGehee made a great leap 
forward in this field in 1974 with his analysis of triple-collision behavior in the 
jj] ■ collinear three-body problem [|9| . He introduced a change of variable to what is now 

known as McGehee coordinates, in which the differential equations for Newton's 
universal law of gravitation are extended to triple collision. By understanding the 
dynamics at triple collision, one understands the dynamics near triple collision via 
continuity of the flow. 

Generating a global analysis of the collinear three-body problem requires an 
understanding of how near-triple collision behavior effects the entire flow. Using 
ideas of Saari and Xia Q , Meyer and Wang undertook this endeavor Q , giving a 
nice picture of the geometry of the phase space for the collinear three-body problem. 

Meyer and Wang's method of analysis involves generating a Poincare slice to the 
flow and partitioning that slice with pieces of stable manifold for triple collision. 
Then in at least some part of the flow the regions generate a subshift of finite type 
on an infinite set of symbols. This is enough to show that chaotic behavior occurs 
in the collinear three-body problem. However, this result raises the question of 

l 
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determining exactly what itineraries of binary collisions are allowed and how the 
set of allowed itineraries changes as the masses vary. 

In order to gain insight on the role of geometry in the collinear three-body 
problem, the author chose to use McGehee coordinates. In contrast, Meyer and 
Wang choose coordinates so that binary collision is represented by a sink. Their 
approach has some advantages; however, in McGehee coordinates binary collisions 
are represented by half-planes. Thus, McGehee coordinates reveal more detail on 
how the stable manifold for triple collision intersects the Poincare slice. 

The crucial construction in this paper is the partition of a Poincare slice into a 
finite number of regions bounded by pieces of stable manifold for triple collision. 
The Theorems in Section 5 give conditions on when this partition is Markov. When 
the partition is Markov, we can construct a semi-conjugacy from the Poincare 
map to a directed graph and exactly describe the set of allowed itineraries. When 
the partition is not Markov, we can construct two sub-shifts of finite type, one 
containing all allowed sequences and one containing guaranteed sequences. These 
two sub-shifts serve as upper and lower bounds on the set of allowed itineraries. 

For details of this approach applied to a model problem in which more detail can 
be explicitly computed, see ||, g, g[ |7| which consider the dynamics of the collinear 
one-bumper two-body problem. 

Except for very special cases, the conditions required by the theorems below can 
only be computed numerically. Even so, the theorems in this paper give a picture 
of what symbolic dynamics occur in the collinear three-body problem and what 
geometric quantities one needs to compute in order to describe those dynamics. 

The symbolic dynamics of the collinear three-body problem are enough to guar- 
antee the existence of special solutions (oscillatory motion) in the ./V-body problem 
(see the Corollary in Section 6). 

2. Hamiltonian Coordinates 

We begin by stating the collinear three-body problem in Hamiltonian coordinates 
and then look at the regularized system in McGehee coordinates. 

Three points masses have masses rrii > and positions q\ < 92 < 93 G R. The 
potential energy is given by 



(i) u = £ 



miirij 



■ 1% - qj 

The motion of the particles under gravitational force is described by the systems 
of differential equations 

(2) m t q t = V qi U, t = 1,2,3 

where \7 qi U is the gradient of U with respect to g, . 

A position, (91,92,93) is called a binary collision if either q± = 92 or 92 =93. If 
91 = 92 = 93 the position is called a triple collision. The above system is defined 
everywhere except at binary and triple collisions. Given an initial positions (not 
at collision) and momenta at time t = 0, a unique solution exists on a maximal 
interval [0,t*). If t* < 00 then the solution is said to have a singularity at t*. 
The only singularities which occur in the collinear three-body problem are due to 
collision fll2|, though non-collision singularities have been found for other TV-body 
problems ]16[| . 
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Double collisions can be regularized. That is, one can change variables so that 
double collision transforms to a regular point of the flow 0] . This extension cor- 
responds to an elastic bounce. McGehee further showed that triple collision is not 
regualrizable. However, in McGehee coordinates, the flow is bounded by an in- 
variant compact manifold which correspond to triple collision. The flow on this 
invariant manifold, called the collision manifold, guides solutions which pass near 
triple collision. 

We change from Hamiltonian to McGehee coordinates via three steps. First 
is a change of variables to polar coordinates. Second is a time scaling change of 
variable so that solutions slow down as they pass near triple collision. Finally, one 
regularizes binary collisions. 

The remainder of this section is a review of the change of variables to McGehee 
coordinates. Readers familiar with this material should feel free to pass on to the 
next section. 

Let q = ((71,(72,92) S K 3 be the vector of positions. Define pi = miqi and let 
p = (pi,P2,P3) & R 3 be the vector of momenta. Let 



M 




Then Equations and @ can be rewritten as 

C/(q) = E 



rrnmj 



.j9t-9i 



(3) Alq = VE/(q). 

We can also write the kinetic energy for the system as 

T(p) = |p T Mp. 
The Hamiltonian for the system is 

H(q,p)=T(p)-C/(q) ) 
and Equation can be written as the system 

q = flp(q,p) = M- X p 

(4) p = -J? q (q,p) = VC/(q). 

The function T is defined everywhere in M 3 . The function U is defined everywhere 
except at collisions. 

Next we break up q into radial and angular components. Define 

r=(q T Afq) 1/2 . 

Notice that a level set of r is an ellipsoid in R 3 . Let S — {q | r — 1}. Then a point 
in S is called a configuration for the system of particles. 
We now define the variables: 

r = (q T Mq) 1/2 

s = r _1 q 

T 
y = p s 

x = p — yMs. 
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In these coordinates, r is the size of the configuration in a inertial norm, s G S is 
the configuration and represents the direction of q, y is the projection of p in the 
direction of q and x represents the direction of change in the configuration. Note 
that s and x are orthogonal. 

In the new polar coordinates, the Hamiltonian, iJ(q, p) = h can be written as 

i(x T M- 1 x + y 2 ) - r^U{s) = h 

and the equations of motion become 

r = y 

y = r" 1 x T M" 1 x-r" 2 f/(s) 

s = r^M^x 

x = -r- 1 yx-r- 1 (x T M- 1 x)+r- 2 U(s)Ms + r- 2 \7U(s). 

We next introduce two new coordinates, u = r 1 ' 2 x and v = r l ' 2 y and scale time 
via dt = r 3 / 2 dt'. In these time-scaled coordinates, the Hamiltonian, iJ(q, p) = h 
can be written as 

|(u T M _1 u + v 2 ) - U(s) = rh 
and the equations of motion become 

r = rv 

v = \v % +u T M- x u-U{s)y 

s = m- 1 vl 

u = -iuu-(u T M _1 u)Ms + l7(s)Ms + VI7(s). 

We reduce the dimension of System ^ by fixing the center of mass at the origin 
and setting the total momentum to zero. This reduces the problem to a four- 
dimensional phase space. 

To express this reduced system, note that there is a unique point on S so that 
Qi = 12 < 13 and the center of mass, Mq = 0. Call this unique point a = 
(«2,&2, a 3)- Likewise, there is a unique point on S so that qi < Q2 = 93 and the 
center of mass, Afq = 0. Call this unique point b = (pi, 62, 62). One can compute 
that < & T Mb < 1. Let A be the least positive solution to 

cos(2A) = a T Mb. 

We now introduce an angular potential function, 

2 



W(s) 



where 



Wi(s) 

W 2 {s) 
Wi(s) 





-(Wi(s) + W 2 (s)+W 3 (s)), 
mim2(l — s) 


sin(2A) 


(b 2 


- bi)Sn(X(l + a)) 

m 2 m 3 (l + s) 




(03 


- a 2 )Sn(X(l - s)) 

Amim3(l — s 2 ) 





and 



(62 - h) sin(A(l + s)) + (a 3 - a 2 ) sin(A(l - s)) 

sin(x) 

Sn{x) — . 
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To regularize double collisions, we introduce a new variable, 

(1 - s 2 )u 



and scale time again via, 



\/W(s) 



dt' = - - dr 



y/W(s) 

Finally, we are ready to write System |] in McGehee coordinates. 
dr A(l - s 2 ) 



dr 



y/W(sj 



dv A 771773^ (1 ' 



Tr = 2^{ l --WW^-±r h) 



ds 



w 



(5) -r- = -sH TTr/ , - fa 2 - 2rfe) + — ^ \ (1 - s 2 - w 2 ) v , ' vw. 

1 ' dr W{s) V ' 2W(s) y ' 2^W{s) 

In McGehee coordinate, the Hamiltonian, H(q, p) = h can be rewritten as 

(6) W{s)(w 2 + s 2 - 1) + (1 - s 2 ) 2 {v 2 - 2rh) = 

System mis defined for all values of r, v, s and w. Note that s is the configuration 
coordinate and varies from —1 to 1. The configuration s = — 1 corresponds to 
91 = 1i < <?3, a left binary collision. The configuration s = 1 corresponds to 
<Zi < <?2 = 93, a right binary collision. Moreover, the vector fields is now defined at 
r = 0, triple collision. 

When r = the Hamiltonian H = h in Equation ra yields the relation 

(7) W(s)(w 2 + s 2 - 1) + (1 - s 2 ) 2 fa 2 ) = 0. 

Equation [7| defines a manifold, M, in R 3 , called the collision manifold, which is 
independent of the total energy, h. The phase space for System || is bounded by 
M. 

3. McGehee Coordinates 

We briefly describe the flow in McGehee coordinates, summarizing results from 
McGehee's research M. For the rest of the paper, we assume that the total energy, 
h, is negative. 

The function, W(s)/(1 — s 2 ) has exactly one critical point at s c . Let 



! W(s c ) 



System H has two equilibria, both of which lay on the collision manifold, M. 
The (r, v, s, w)-coordinates of these two equilibria are c = (0, — u c ,s c ,0) and d = 
(0, v Cl sc, 0) (see Figure |l|) 

On the collision manifold, c and d each have a one-dimensional stable manifold 
and a one-dimensional unstable manifold. The flow on the collision manifold is 
gradient-like with respect to level sets of v (that is, dv/dr > on M). A sketch of 
these stable and unstable manifolds is given in Figure |[ 
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Figure 1. Collision manifold with stable and unstable branches 
of d on M. 





Figure 2. Ejection-collision orbit with stable and unstable 
branches of c and donM. 



There is a heteroclinic connection between c and d corresponding to an ejection- 
collision orbit (see Figure g). This solution is homographic, that is, the configura- 
tion is constant. So, along the ejection-collision orbit, the s-coordinate is constant, 
s c . The forward solution of the orbit limits onto c and the backwards limit is d. 
Thus the ejection-collision orbit begins and ends at triple collision and does not 
pass though a binary collision. 

We say any solution which has a forward limit on c, "ends in triple collision" 
and we say any solution which has a backwards limit from d, "begins in ejection" . 
The set of solutions which end in triple collision forms the two-dimensional stable 
manifold of c, W s (c) and the set of solutions which begin in ejection forms the 
two-dimensional unstable manifold of d, W u (d). Any orbit in the intersection of 
W s (c) and W u (d) is called an ejection-collision solution. 
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Binary collisions occur at s — ±1. The energy relation in Equation p] requires 
w = at binary collision. Thus we must use both the r and v coordinates distinguish 
one binary collision from another. Although the System appears to separate the 
r-coordinate from the others (and indeed does so for \s\ < 1), we will have to include 
the r-coordinate at binary collisions. 

All solutions, except for the nomographic ejection-collision orbit, pass through 
at least one binary collision. This makes the pair of half planes, 

r = {(r, v, s, w) | r > 0, \s\ — 1, w = 0} 

the appropriate choice for a Poincare slice. Thanks to work by Mayer and Wang, 
we know something about the geometry of the intersection of the stable manifold 
for triple collision and this Poincare slice, T. 

The motivation for studying how the stable manifold intersect the Poincare slice 
is given by the following argument. Any point on the Poincare slice has a forward 
itinerary of binary collisions, perhaps ending in triple collision. Given a continuous 
arc on the Poincare slice whose endpoints have different itineraries, there must be 
some point on the arc which is also on W s (c). That is, W s (c) divides the Poincare 
slice into regions with different itineraries. 

We distinguish the two half-planes of the Poincare slice by L for s = — 1, and R 
for s — \. The intersections of W s (c) with L and R are arcs with two endpoints 
on r = or loops with one endpoint on r = 0. We label these arcs and loops by 
their itineraries using L's and R's for binary collisions and C for triple collision. We 
use the subscript + to designate the location of the set with a given itinerary. For 
example the itinerary LL+RC is the set of initial conditions on L which in forward 
time pass through R once before triple collision and in backwards time pass through 
L again. The itinerary LLR+C, however, designates the set of initial conditions on 
R which lead directly to triple collision and whose prior two binary collisions were 
on L. 

Solutions on M, may lay on the stable manifold for c or for d. Such solutions will 
terminate with the symbol c or d as appropriate. For example, L+Rc designates the 
initial condition on M which traverses to R and then limits onto c without passing 
though another binary collision. Likewise, L+Ld designates the point on M which 
begins on L, returns to L and then limits onto d without passing though another 
binary collision. Using this itinerary notation for the intersections of W s (c) and 
the Poincare slice, we next summarize Meyer and Wang's results |ll| in the setting 
of McGehee coordinates. 



There is a unique arc on L with itinerary L+C (see Lemma 4.1). This arc has 
two endpoints on M, one with the itinerary L+c and the other L+d. Likewise, 
there is a unique arc on R with itinerary R+C. Its endpoints have itineraries R*c 
and R±d (see Figure ||). We know these arcs are unique since they must each be 
homotopic to the ejection-collision orbit between c and d which passes through no 
binary collisions. 

To designate a solution which begins in ejection we use the symbol E. For 
example, ELR+LC designates the initial condition(s) on R which pass though L 
and limit onto c. These initial condition(s) must also in backwards time pass 
through L then limit onto d without passing through any other binary collisions. 

Intersections of the ejection manifold, W u (d), with the Poincare slice have a 
simple relation to intersections of W s (c) because System || is reversible. That is, it 
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Figure 3. First pullback of stable manifold for triple collision. 



has a time symmetry. Let 



•</• 



-w 



If (r,v, s, w) is a solution in time r, then so is (r, v, S, w) in time r. 

This means that the arcs L*C and EL* are reflections of one another over the 
v = axis. In general the arc or loop with itinerary ao.a1a2--.akC flips about the 
V = axis to an arc or loop with the itinerary Eak---a2aiaQ. and vice versa. 

With the itinerary notation, we can briefly describe the work of Meyer and 
Wang. They generate their main results by focusing on regions whose itinerary 
begins L+L, L+RL, L+RRL, and in general L*R n L (where R n denotes a string of 
n R's) . These regions are bounded by pieces of stable manifold whose itinerary is 
LR n C. They show that the first return of these regions to L generates a sub-shift 
on an infinite number of symbols. Since not all regions are guaranteed to intersect 
the first returns, the exact dynamics can not established although it is clear that 
the dynamics are rich. 



4. The Stable Manifold of Triple Collision 

In this section we characterize pullbacks of the stable manifold for triple collision, 
W s (c), on L and R. Recall we are studying the dynamics of the System ||. We 
have chosen a Poincare slice, T, which is transverse to the flow || made up of two 
half-planes, L and R, corresponding to left and right binary collisions. The flow 
induces a map on T. We denote this Poincare (first return) map by P. We now 
want to show how the stable manifold for triple collision, W s (c), intersects T. 
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Lemma 4.4 




Lemma 4.2 



Lemma 4.3 



Figure 4. Typical arcs on L labeled with corresponding Lemmas. 

Lemma 4.1. The first pull backs of the collision manifold, L+C and R+C , are 
smooth arcs with end points L+c, L*d and R*c, R+d. 

Proof. Since the flow is smooth, L+C and R+C are deformations of the ejection- 
collision orbit between c and d. Hence L+C and R+C are smooth arc with end 
points on M. Since the endpoints must limit to c and d, the endpoints of L*C are 
L+C and L+d. Likewise, the endpoints of R+C are i?*c and R+d. □ 

We now show how to continue pulling back pieces of stable manifold for triple 
collision (see Figure 0). The following Lemmas are consistent with Meyer and 
Wang's results. However, Meyer and Wang proved them for the case where left 
binary collisions were not regularized. Here, we show them for the fully regularized 
system. 

Lemma 4.2. // an arc, 7 in W s {c) D T has endpoints x and y on M and the arc 
does not intersect EL* or ER* then the pullback of"/, is an arc on T with endpoints 
P- I (x) and P- I {y) onM. 

Proof. Clearly, x pulls back to P^ 1 (x). Points on 7 near x pull back near P^ 1 (x). 
The only obstruction to pulling back a point on 7 to T is if the point pulls back to 
ejection. Since 7 has no such obstruction, all of 7 pulls back to T as described. 

□ 



Lemma 4.3. If an arc in W s (c) n T has endpoints x and y on M and the arc 
intersects EL* or ER+ then the segment, 8, from x (or y) to the first intersection 
with EL* or ER* pulls back to an arc on V with endpoints P~ l (x) (or P~ l (y)) 
and either L*d if the segment pulls back to L or R*d if the segment pulls back to 
R. 

Proof. Assume that S has an endpoint at x. Clearly, x pulls back to P~ 1 (x). 
Points on 8 near x pull back near P~ 1 (x). We continue pulling back points along 8. 
Points near the intersection of 8 and EL* or ER+ must pull back arbitrarily close 
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L 



R 



EL J 




Figure 5. An arc and its pullback for Lemma 4.2 



EL* 




b> v 



P- 1 ^)! I EL, 




P- I (y) 



R 

R*d 



p-HSi) 



P- I {x) 



Figure 6. Two arc segments and their pullbacks for Lemma 4.3 



to the equilibrium, d, hence return to T arbitrarily close to one branch of the stable 
manifold for d on M. Thus points on 8 near the intersection of 8 and EL* or ER* 
must pull back to T arbitrarily close to either L±d or R±d. By continuity, points 
near the intersection of 8 and EL+ or ER+ must pull back to the same half-plane, 
L or R, as the rest of S. 

□ 



Lemma 4.4. If a segment of an arc, 8, in W s (c) n T has both endpoints on EL+ 
or ER± and has no other intersections with EL+ or ER+, then the segment pulls 
back to a loop on T with both endpoints at either L+d if the segment pulls back to 
L or R±d if the segment pulls back to R. 

Proof. Since 8 intersects EL+ or ER+ only at its endpoints, all of 8 pulls back to 
either L or R. Points on 8 near the intersections of 8 and EL+ or ER+ must pull 
back arbitrarily close to the equilibrium, d, hence return to T arbitrarily close to one 
branch of the stable manifold for donM. Thus points on 8 near the intersections 
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EL*\ 




L 



L*d R+d 



R 



Figure 7. Two arc segments and their pullbacks for Lemma 4.4 



of 6 and EL* or ER+ must pull back to T arbitrarily close to either L+d if the rest 
of <5 pulls back to L or R+d if the rest of <5 pulls back to R. 

Pullbacks of loops are covered in earlier Lemmas since we may let x = y in 
Lemma 4.2. □ 



By the lemmas above, we know that the endpoints of arcs and loops in W s (c) PlT 
are on either the stable manifold of c or d on M. By McGehee, we know that the 
stable manifold for c has two branches on M. The stable manifold for d has two 
branches. In backwards time, one branch passes first through L and then alternates 
between R and L, eventually going down one of the legs. The other branch of W s (d) 
on M in backwards time passes first through R and then alternates between L and 
R, eventually going down one of the legs (see Figure [l]). 

One branch of W s (d) on M has an itinerary ...LLLRLRL...Ld. We denote the 
length of the alternating part of this itinerary £m- Likewise, we denote the length 
of the alternating part of the branch of W s (d) ending ...Rd by tm- For example, a 
branch of W s (d) on M with itinerary ...LLLRLRLd yields £m. = 5. 

Since the branches of W s (d) can not intersect we have the relation 



(8) 



Km -nyr| < 1- 



If £m. = ?"m then either each branch goes down a different leg of M or there is a 
heteroclinic connections between c and d on M, i.e. W s (d) = W u (c). Otherwise 
£m and tm differ by one and both branches must go down the same leg. 

McGehee shows that the v-coordinate along flow on M is non-decreasing. Meyer 
and Wang show that the stable branches for d on M must end with RLd or LRd. 
That is, 

t-w, ^m > 2. 
These two facts have the following geometric consequence. 

Lemma 4.5. The stable manifold of c with itinerary L+C and the unstable mani- 
fold with itinerary EL+ intersect on L. Likewise, R+C and ER+ intersect on R. 

Proof. We only look at the case of the intersection of L*C and EL* on L since the 
other case is similar. 

Let v{x) be the ^-coordinate of a point on M. Since v is non-decreasing along 
solutions on M, we know from the definition of stable and unstable manifolds that 
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the following inequality holds 

v{L±C) < v(c) < < v{d) < v(dL*) 

If v(L+d) < v{cL+), then the inverse image of L+d under the Poincare map must 
be on L since the pullback of L+d can not intersect the trajectory between c and 
cL*. This means that one branch of the stable manifold for donM would have the 
itinerary LLd which contradicts Meyer and Wang's result. Thus v(cL+) < v(L+d). 
Since v{cL+) + v(_L*oQ = we have the inequality 

(9) u(£*C) < v{c) < v{cL+) < < v(L*cC) < v(d) < v{dL+) 

The continuous arc L+C has endpoints at L+c and L+d. The continuous arc EL+ 
has endpoints at dL+ and cL+. By the inequality in Equation |9j, L+C and EL+ 
must intersect at least once. □ 



The segment of L±C from L±d to the first intersection of EL* by Lemma 4.3 
must pullback under P^ 1 to an arc on R with endpoints at R+Ld and R*d. If this 
arc does not intersect ER+ then we continue pulling back this segment until there is 
a first intersection with either EL+ or ER+ . Denote by i n the number of pullbacks 
required for the first intersection. Likewise, let r n denote the number of pullbacks 
required for the segment of R+C from R+d to the first intersection with ER+ to 
next intersect EL+ or ER+. 

5. Main Results 

The Lemmas of the previous section establish that W s (c) n L is made up smooth 
arcs and loops whose endpoints are on the stable manifolds of c and d on M. Thus 
W s (c) (~1 r generates an infinite number of regions on P. However, only a finite 
number of regions are needed to determine or bound the symbolic dynamics. 

The values of £m., Tm., £n, and r n characterize the global dynamics. The following 
theorems connect £m , fm , in , and r n with the global dynamics. 

Before proceeding, we review the definition of a sofic system. A finite directed 
graph whose arrows are labeled and the labels may be used more than once is called 
a sofic system. The dual to a sofic system, then, is a finite directed graph whose 
vertices are labeled and whose labels may be used more than once. Sofic systems 
are a generalization of sub-shifts of finite type. See |15[ and M for more information 
on sofic systems and their role in dynamical systems. 

Theorem 5.1. If for three given masses £m. = ?*M = ^n + 1 = r n + 1 then the 
set of allowed itineraries for the system is given by the dual of a sofic system (see 
Figures \1(\ and \li 



Proof. The arc EL+ divides L into two regions. We call the bounded region the 
inside of EL+ and the unbounded region the outside of EL+. Likewise, we define 
the inside of ER± and the outside of ER+ 

We next divide the inside of EL+ into a finite number of regions via the following 
procedure. The segment(s) of L+C inside EL+ is (are) the primary left segment(s), 
denoted PLS. Note that the segment of L+C from from L±d to the first intersection 
with EL+ is always contained in the the PLS. If PLS contains any other segments, 



they must be of the form described by Lemma 4.4. 

By the definition of in, the first in — 1 pullbacks of PLS are (or at least contain) 
arcs inside ER+ and EL*. We call these arcs PLS~ l ,PLS~ 2 ..., PLS~ in+1 succes- 
sively. Likewise, we define PRS and its pullbacks, PRS- 1 , Pi?5" 2 ..., PRS- ln+l . 
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PLS^ 



PRS 



L+d 



P- X {PRS)[ 



EL*\ P- 2 (PLS)[ 



R+d 



)P- 1 {PLS) 
P-^L+d) 

p-\R+d) )p-i(PRS) 



\ER* 





P- 



i+i 



{PRS) 



p-*n+l(pLS) 



FIGURE 8. Dividing L and R each into £ n 
labeled for £ n even. 



2 regions. Diagram 



The arcs, EL± : ER* : PLS and its pullbacks and PRS and its pullbacks divide L and 
R each into £ n + 2 regions. Figure |8| describes the case where PLS and PRS each 
only have one segment. If they had additional segments, loops would be attached 
at the pullbacks of L+d and/or R+d. 

Points on M of the form p- 2l (L+d) or p- 2i+1 (R*d) (i > 0) which are in the 
closure of the inside of EL+ are called inside points of EL+. Likewise, points on 
M of the form p- 2l (R±d) or p- 2i+1 (L±d) {i > 0) which are in the closure of the 
inside of ER+ are called inside points of ER+. Points on M of the form P~ l (R+d) 
or P^^^L+d) which are not inside points are called outside points. 

We define a vertical strip as a simply connected subset of a region whose inter- 
section with the boundary of the region is made up of one of the following: 

1. at least one inside point and a segment of EL+ or ER-,, 

2. at least one outside point and a segment of EL+ or ER+ 

3. at least two inside points 

4. at least two outside points 

The pullback of an inside point is either an inside point or an outside point. The 
pullback of an outside point is an outside point. The pullback of a region arbitrarily 
near EL+ or ER+ is a region either arbitrarily near L+d or R+d. So the pull back 
of a vertical strip at least contains vertical strips. If the pullback of a vertical strip 
crosses over EL+ or ER+, two segments of the pullback are vertical strips while the 
rest of the segments only have ends on EL+ or _Ei?*, and hence are not vertical (see 
Figure ||) . 
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Figure 9. Shaded vertical strip on L and its pullback on R. 



Although we can make a directed graph showing how regions map to one another 
under -P -1 , we must also show that the pull back of a region's vertical strip contain 
vertical strips in the pullback of the region. 

We next describe the pullbacks of various cases of vertical strips. Our goal is to 
show that the division of L and R described above is a Markov partition of T for 
P^ 1 (hence Markov for P by reversability) . 

We begin with a vertical strip inside EL* and above the arc of PLS from L*d 
to the first intersection with EL*. This vertical strip pulls back to a vertical strip 
between R*d and RL*d inside PLS -1 . Further pulls backs generate vertical strips 
in regions bounded PLS~ % for 1 < i < l(~\ ~ 1- 

A vertical strip inside PLS ln+1 contains at least two vertical strips. One strip 
includes a segment of EL* or ER* and the outside point, P~ tn (L*d). The other 
strip includes EL and the inside point P~ ln+1 (R*d). This second strip, though 
inside EL* or ER* is also outside the regions PLS~ % or PRS~ % for 1 < % < l^ — 1. 

A vertical strip from EL* or ER* to any outside point has a pullback which 
includes a vertical strip from inside EL* or ER* and above either the arc of PLS 
from L*d to the first intersection with EL* or the arc of PRS from R*d to the first 
intersection with ER* . 

Like a vertical strip inside PLS £n+1 , a vertical strip from a segment of EL* or 
ER* and with P~ £n+1 (L*d or P~ tn+1 {R*d) contains at least two vertical strips. 
One strip includes a segment of EL* or ER* and the outside point, P~ £n (L*d). 
The other strip includes EL and the inside point P~ in+1 (R*d). This second strip, 
though inside EL* or ER* is also outside the regions PLS~ l or PRS~ l for 1 < i < 
in - 1. 
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Figure 10. Directed graph for Theorem 5.1 where £ n is even. 



The above argument can be repeated beginning with a vertical strip inside ER+ 
and above the arc of PRS from R±d to the first intersection with ER+. 

Since vertical strips in a region pull back to vertical strips in the pullback of the 
region, the partition of L and R given above is Markov. If we label the regions, then 
there is a subshift of finite type which semi-conjugate to P^ 1 on T. By reversability, 
we can switch the direction of the arrows and get a subshift of finite type which is 
semi-conjugate to P on I\ 

Replacing region names with "L" for regions on L and "i?" for regions on R 
yields a directed graph on two symbols (the dual of a sofic system). This directed 
graph in general can be simplified. The simplified graph describes all itineraries of 
binary collisions. 

To see this last claim, any initial condition on L or R is in some region as 
defined above. Thus the point maps about the regions according to some path in 
the sub-shift, so its itinerary is described by the graph. 

Conversely, given any itinerary described by the graph, there is at least one path 
in the sub-shift which accomplishes the given itinerary. Choose a closed vertical 
strip in the first region for such a path. Since the pullback of vertical strips includes 
contains vertical strips we can generate a nested sequence of closed vertical strips 
which obtain arbitrarily many terms in the desired sequence. The infinite intersec- 
tion of these closed and nested vertical strips is non-empty, thus guaranteeing at 
least one point which achieves the given itinerary. □ 



Carrying out the procedure outlined in Theorem 5T we see that for the case 
£m = ?"m = ^n + 1 = fn + 1; the generated dual to a sofic system is given by 
Figures pi] and O. Note that the symbol A£ means an alternating sequence of 
I/'s and i?'s beginning with L of length £ n — 1. Likewise, the symbol A£~ x means 
an alternating sequence of L's and -R's beginning with R of length £ n — 1- 

Note that for ^m = ?'m = Cn + 1 = ^n + 1 = 2 the dual to a sofic system reduces 
to the full shift on two symbols. 

Theorem 5.2. If for three given masses £m = tm + 1 = ^n + 2 < £ n or tm = 

C-tsa + 1 = C-n + 2 < r n then the set of itineraries for the system is bounded between 
two duals of sofic systems (see Figures uj to UJa). 



Proof. We will assume the first case, £m 
for the second case is similar. 



?"m + 1 = ^n + 2 < in i since the argument 



Using the notation from Theorem 5.1 we divide the inside of EL+ and ER* 



into a finite number of regions bounded by PLS % for < i < r n and PRS l for 
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Figure 11. Directed graph for Theorem 5.1 where £ n is odd. 



< i < rp, — 1. We also make regions outside ER+ or EL* bounded by PLS l for 
r n + 1 < i < in ~ 1- 



By the argument in Theorem 5.1, a region's vertical strip pulls back to vertical 
strips in the region's pullbacks with one exception. A vertical strip inside PLS~ in+1 
will not pull back to a vertical strip inside EL* or ER* even though the region 
bounded by PLS~ in+1 and M must cross either EL* or ER*. This means that a 
directed graph showing how regions map under P" 1 will not be Markov. However 
a point on T must move from region to region according to the the directed graph 
so the directed graph contains all allowed itineraries of regions even though some 
itineraries may not be achieved. 

So we begin with the directed graph showing how regions map under P , By 
reversability, switching the direction of the arrows we have a directed graph showing 
how regions map under P. Replacing region names with "L" for regions on L and 
"i?" for regions on R yields a directed graph on two symbols (the dual of a sofic 
system). This directed graph in general can be simplified. The simplified graph 
contains all itineraries of binary collisions. 

To bound the set of allowed itineraries from below, we assume that pullbacks 
of PLS never intersect EL* or ER*. That is we assume r n = +oo. Although 
this is not the case, the condition guarantees a region's vertical strips pull back to 
vertical strips in the pullback of the region. That is, the directed graph showing 
how regions map under P~ x with the condition that r n = +oo is Markov, hence 
the sub-shift contains orbits which are guaranteed to occur in the actual system. 

We begin with the directed graph showing regions map under P^ 1 with the 
condition that r n — +oo. By reversability, switching the direction of the arrows 
we have a directed graph showing how regions map under P. Replacing region 
names with "L" for regions on L and "i?" for regions on R yields a directed graph 
on two symbols (the dual of a sofic system). This directed graph in general can 
be simplified. The simplified graph contains all guaranteed itineraries of binary 
collisions. □ 



Carrying out the procedure outlined in Theorem 5.2 we see that for the case 
£m = ?"m + 1 = ?"n + 2 < In, the generated duals to sofic systems is given by 
Figures [l2| to 15. For the case tm = ^m + 1 = ^n + 2 < r n , the generated duals to 
sofic systems is given by Figures O to n5^ after exchanging L's and i?'s. 



Theorem 5.3. If for three given masses, if £m. , ^m , ^n and r n do not meet the 
criteria of Theorem 5.1 or Theorem 5.i then the set of itineraries for the system is 



bounded between two duals of sofic systems. 
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FIGURE 12. Directed graph for Theorem |5.2| containing allowed 
sequences for £-m even. 




Figure 13. Directed graph for Theorem p.2| containing allowed 
sequences for £-m odd. 



Proof. If ^m> r M, ^n and r n do not meet the criteria of Theorem 5.1 or Theorem 5.2 



then (using the notation of Theorem 5.1) one of the pull backs of PLS or PRS 
must cross EL* or ER± when both of cndpoints of that arc are either inside EL* 
or ER+. This means that a vertical strip which landed in this region will pass to a 
non- vertical strip outside EL* or ER±, even though the region pulls back outside 
EL* or ER+. This means a Markov partition is not possible f or a directed graph 
on any finite set of regions similar to Theorem ^J or T heor em 5.2 . 

If I'm = Hvi then the directed graph from Theorem 5T for these values of ^m 
and rjvi must contain guaranteed dynamics for our system since the premature 
intersection described above only adds to the possible set of itineraries. 



Likewise, if ^m ^ tm. then the associated directed graph from Theorem 5.2 which 
contains guaranteed sequences also contains guaranteed sequences for our system 
since the premature intersection described above only adds to the possible set of 
itineraries. 

To get an upper bound for our system we divide L and R into a finite number 
of regions bounded by EL+, ER+, PLS- 1 for < i < £ n - 1 and PRS- i for 
< i < r n — 1. Each of these pullbacks are single arcs by the conditions on £ n and 
7'n- The directed graph describing how these regions map under P _1 must contain 
all allowed itineraries of regions even though some itineraries may not be achieved. 
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Figure 14. Directed graph for Theorem 5.2 containing guaran- 
teed sequences for ^m even. 




R <- 



R 





Figure 15. Directed graph for Theorem 5.2 containing guaran- 
teed sequences for £m odd. 



So we begin with the directed graph showing how regions map under P . By 
reversability, switching the direction of the arrows we have a directed graph showing 
how regions map under P. Replacing region names with "L" for regions on L and 
"R" for regions on R yields a directed graph on two symbols (the dual of a sofic 
system). This directed graph in general can be simplified. The simplified graph 
contains all itineraries of binary collisions. □ 



Carrying out the procedure outlined in Theorem 5J3 as an example, we see that 
for the case £m. = 5, tm = 6, £ n — 2 and r n = 7, the generated dual to a sofic 
system containing all allowed sequences is given by Figure |16| . 

Note that for £ n = r n = 1 the upper bound for the dynamics is the full shift on 
two symbols, hence a trivial upper bound. 



6. Oscillatory Motion 

One application of the Theorems in Section 5 is to the presence of oscillatory 
motion in the iV-body problem. Saari and Xia explored possible behaviors in the 
TV-body problem as t — > oo |L4j. One such behavior is oscillatory motion, that is, 
a mutual distance coordinate r, so that as t — > oo, the limsupr = oo while the 
liminf r < oo. 

For the collinear three-body problem, oscillatory motion requires that either 
mi or TO3 takes longer and longer excursions for the other binary pair, each time 
returning to the binary pair before its next excursion. Saari and Xia showed that 
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Figure 16. Directed graph for Theorem ^_^ containing all allowed 
sequences for the case £m = 5, tm = 6, £ n = 2 and r n = 7. 



it is enough to prove the presence of itineraries of the form 



.L a \..L a 



.R ai ...R a2 ...R a3 ...R ai . 



so that a, — > oo as i — > oo. The existence of such itineraries guarantees the existence 
of oscillatory motion. 

Saari and Xia study the return map on the zero momentum set for m^. To 
guarantee a return map, they restrict the masses to the case where triple collision 
can lead to arbitrarily high velocities. This is exactly when the stable manifold of 
d and and the unstable manifold of c intersect transversely. In this case, the set if 
allowed itineraries is given by the full-shift on two symbols (same as £~m = tm = 

Once they establish the existence of oscillatory motion in the collinear three- 
body problem for some sets of masses, they note that the motion can be extended 
to the iV-body problem. They conclude that for any N > 3 there exist masses and 
initial conditions so that oscillatory motion exists. 

From the Theorems in Section 5 it is clear that for all sets of masses which admit 
transverse intersections of the stable manifold of d and and the unstable manifold 
of c, itineraries of the form 

...L ai ...L a2 ...L a3 ...L a4 ... or ...R ai ...R a2 ...R a3 ...R ai ... 

so that at — > oo as i — > oo exist. We are thus led to the same conclusion as a 
Corollary to our Theorems. 

Corollary 6.1. For the N-body problem, N > 3, there exist positive masses and 
initial conditions so that oscillatory motion occurs. 
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